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SUMMARY 

Understanding viral transmission dynamics within populations of reservoir hosts can facilitate 
greater knowledge of the spillover of emerging infectious diseases. While bat-borne viruses are of 
concern to public health, investigations into their dynamics have been limited by a lack of 
longitudinal data from individual bats. Here, we examine capture-mark-recapture (CMR) data 
from a species of Australian bat ( Myotis macropus) infected with a putative novel 
Alphacoronavirus within a Bayesian framework. Then, we developed epidemic models to estimate 
the effect of persistently infectious individuals (which shed viruses for extensive periods) on the 
probability of viral maintenance within the study population. We found that the CMR data 
analysis supported grouping of infectious bats into persistently and transiently infectious bats. 
Maintenance of coronavirus within the study population was more likely in an epidemic model 
that included both persistently and transiently infectious bats, compared with the epidemic model 
with non-grouping of bats. These findings, using rare CMR data from longitudinal samples of 
individual bats, increase our understanding of transmission dynamics of bat viral infectious 
diseases. 

Key words: Bayesian analysis, coronavirus, epidemiological modelling, Myotis macropus , persistent 
infection. 


INTRODUCTION 

Coronaviruses have been increasingly recognised as a 
human public health issue following the emergence of 
high-impact zoonotic diseases from bats - the mam¬ 
malian order that hosts the largest diversity of 
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coronaviruses [1], Examples include the coronaviruses 
that caused severe acute respiratory syndrome 
(SARS), which claimed 916 lives out of 8422 cases 
from November 2002 to August 2003 [2], and 
Middle East respiratory syndrome (MERS), which 
has claimed 608 lives out of 1449 cases as of 31 
August 2016 [3]. SARS-like coronaviruses (SL-CoV) 
are maintained in bats [4], and MERS coronavirus 
(MERS-CoV) is assumed to have originated in bats 
[5]. The findings that genetically diverse SL-CoV 
strains share high similarity with SARS coronavirus 
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(SARS-CoV) and that bats harbour diverse corona- 
viruses which can be classified to the same coronavirus 
species as MERS-CoV suggest that recurrent corona- 
virus epidemics and pandemics in humans are likely 
[6]. Genetic analyses of coronaviruses in bats sug¬ 
gested that the diversity of coronaviruses in bats 
may provide further opportunities for spillover into 
other species [7]. Mitigation of spillover of the viruses 
is based on understanding of the maintenance mech¬ 
anism of multihost viruses in reservoir hosts. 

Bats may have idiosyncratic immune responses, allow¬ 
ing them to be infected (with no signs of disease) with 
viruses that are highly pathogenic in other species [8], 
This feature of bats may contribute to their ability to 
host zoonoses, including SARS, MERS, Ebola, 
Marburg, Nipah and Hendra viruses. A number of 
hypotheses have been proposed to explain this: 
some authors suggest that bats' ability to fly induces 
fever, helping to control viruses [9]; others suggest that 
immune system adaptation allows tolerance of intracellu¬ 
lar pathogens [10] and that co-evolution during a long 
history between specific viruses and bat hosts results in 
no pathogenicity to bats whereas recent host shifts from 
bats to other mammals result in high fatality rates [11,12], 

Two conditions must be met for bats to be reservoir 
hosts of emerging infectious diseases: first, the ability to 
maintain pathogens in their populations, and second, 
the ability to transmit those pathogens to another 
species. This paper concentrates on the first of these con¬ 
ditions. Features that may contribute to viral mainten¬ 
ance in bat populations include: (1) a metapopulation 
structure that avoids viral extinction across the total 
population via sub-populations, which allow reinfection 
through movement [1 3,14]; (2) reduced metabolism dur¬ 
ing hibernation that facilitates viral survival [1 5]; and (3) 
persistent infections with chronic shedding or intermit¬ 
tent recrudescence [16, 17]. This study focuses on the fea¬ 
ture of persistent infection. Persistent infections of 
SARS-CoV in vitro studies support the plausibility of 
persistent infections of coronaviruses in bats [18, 19]. In 
persistent infections, a virus is not cleared from the 
host but remains associated with specific cells for a 
long period. Regardless of the mechanisms operating 
within hosts, persistent infection contributes to viral 
maintenance within populations and has been consid¬ 
ered as a mechanism of viral maintenance in reservoir 
hosts of other emerging infectious diseases [20]. 

This study is based on data obtained in a previously 
reported capture-mark-recapture (CMR) study [21] 
of a maternal roost of Myotis macropus. This microbat 
(Microchiroptera ), also known as large-footed myotis, 


Table 1. The CMR data composition of coronavirus in 
52 M. macropus 

Recapture 


Coronavirus RNA detection 

Yes 

No 

Total 

Multiple 


0 

7 

Single 

16* 

5 

21 

No 

19 

5 

24 

Total 

42 

10 

52 


Seven ‘persistently infectious bats’ were recaptured bats with 
multiple detections of coronavirus RNA, and 16 ‘transiently 
infectious bats’ were recaptured bats with a single detection 
of coronavirus RNA. Twenty-three ‘infectious bats’ were 
persistently or transiently infectious bats. ‘Infectious bats’ 
did not include five bats, which were not recaptured, with 
a single detection of coronavirus RNA. 

* Seven recaptured bats with multiple detections of corona¬ 
virus RNA were identified with total 18 detections of 
coronavirus RNA. 
f Persistently infectious bats. 

J Transiently infectious bats. 

is widely distributed in Australia [22], M. macropus 
inhabits areas close to waterways in small groups [21], 
foraging on aquatic invertebrate and small fish [22]. 
The bats in the study region fonn maternity roosts bian- 
nually in October and January, and deliver single young 
[21]. Gestation and lactation periods are 12 and 8 
weeks, respectively [21]. 

Smith undertook a CMR study of M. macropus in 
the lifting holes of a bridge in southeast Queensland, 
Australia (Supplementary Table SI) [21], Full details 
of the methodology are described in the original 
study [21], Briefly, the CMR data of 52 Australian 
bats ( M. macropus ) were collected during nine captur¬ 
ing occasions over 12 weeks from 13 January to 31 
March 2009, which overlapped with the January 
breeding season and subsequent lactation period for 
this species [21]. Coronavirus RNA was detected in 
faecal pellets or anal swabs using RT-PCR (reverse 
transcription polymerase chain reaction) targeting a 
conserved region of a coronavirus gene [21]. Of the 
52 unique individuals, 42 were recaptured at least once, 
and of these, 7, 16 and 19 bats respectively showed mul¬ 
tiple detections, a single detection, and no detection of 
coronavirus RNA (Table 1). The data from bats with 
multiple detections of coronavirus RNA were suggestive 
of persistent infections [21], The CMR dataset are 
unusual and particularly valuable in that the data con¬ 
tain individual tracking records with infection states. 
Longitudinal sampling of individual bats is necessary 
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Fig. 1. CMR data analyses across eight recapturing occasions, (a) Survival and recapture rates. Black and white circles 
represent survival and recapture rates, respectively, (b) Transition rates between uninfectious and infectious states. Black 
and white circles represent transition from uninfectious to infectious state and from infectious to uninfectious state, 
respectively. Error bars indicate 95% CrI. 


to explore persistent infection and to test whether persist¬ 
ent infection is a possible viral maintenance mechanism 
in a bat population. For many bats, large population size 
and migratory behaviour have impeded obtaining indi¬ 
vidual data. Most field studies have therefore been cross- 
sectional [14, 23, 24] rather than following individuals 
through time, limiting development of understanding 
the maintenance mechanism of bat-borne viruses. 

Here, we apply quantitative analyses in a Bayesian 
framework to the CMR data from Smith [21] to derive 
parameter estimates for epidemic models to investi¬ 
gate the effects of persistent infection on viral main¬ 
tenance in a bat population. 

METHODS 
CMR data analyses 

We analysed the CMR data using multistate models to 
estimate survival rates (<p), recapture rates (p), and tran¬ 
sition rates (i//) (Fig. 1). The multistate models included 
two states, infectious and non-infectious, with transitions 
in both directions [25]. Each time step was a week. In 
addition, to investigate persistent infection in the bat 
population, we explored multistate models in which 
bats were divided into two groups, based on the fre¬ 
quency of coronavirus RNA detection. Recaptured 
bats with a single detection were referred to as transiently 
infectious bats and recaptured bats with multiple detec¬ 
tions were referred to as persistently infectious bats. 
Bats that were only captured once were excluded from 
analysis (Table 1). Five bats were once captured with a 
single detection of coronavirus RNA, and those five 
bats were excluded from two groups. 


To analyse the CMR data, we chose a Bayesian 
method over a frequentist approach because the 
Bayesian method is more appropriate in dealing 
with a small amount of data, relying less on large sam¬ 
ple asymptotic approximations [26], The Bayesian ana¬ 
lyses were conducted in R [27] and in OpenBUGS using 
the R package ‘R20penBUGS’ [28]. OpenBUGS was 
used to run three independent chains of an MCMC 
(Markov chain Monte Carlo) sampler for 10 000 itera¬ 
tions each, after discarding the initial 1000 samples as a 
‘burn in’. The mean of each of the parameters was cal¬ 
culated, as were the 2-5th and 97-5th percentiles of the 
parameter distributions (95% Bayesian credible inter¬ 
vals (CrI)). The relevant R code is provided as 
Supplementary material (R code SI and S2). 

Multistate model selection 

We used CMR multistate model selection to determine 
whether ‘grouping’ of persistently and transiently infec¬ 
tious bats and the inclusion of infectious and non- 
infectious states were supported by the CMR data. 
The grouping and the multiple states were applied to 
survival, recapture and transition rates of the CMR 
data. Comparisons between the candidate multistate 
models were assessed with deviance information criter¬ 
ion (DIC) [29], and the most parsimonious model was 
selected for subsequent simulated epidemic models. 

Parameterisation 

Survival rates ( cp ) were used to calculate mortality 
rates (/<). Transition rates (<//) from the infectious to 
non-infectious state (y/ ui ) were used to calculate 
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Table 2. Parameters of coronavirus infection in M. macropus; each model time step is 1 week 


Parameter 

Infection 

type 

Symbol 

Estimate or range 

Source 

Transmission rate 


P 

0-0104 

Calculated from equations for basic 





reproductive rate [31] 

Recovery rate 

All 

y 

0-5638 (95% CrI 

Transition rates from infectious state to 




0-3236-0-8031) 

non-infectious state (y/ ui ) [21] 


Persistent 

)’p 

0-3354 (95% CrI 





0-1210 0-6518) 



Transient 

7t 

0-8582 (95% CrI 





0-4985-0-9943) 


Rate of waning immunity 


CO 

0-0833 

[30] 

Mortality rate in 

All 

/^u 

0-0152 (95% CrI 

1-survival rate [21] 

non-infectious state 



0-0732-0) 



Persistent 

j^up 

0-0617 (95% CrI 





0-2817-0-0017) 



Transient 

/^ut 

0-0166 (95% CrI 





0-0815-0) 


Mortality rate in infectious 

All 

Mi 

0-0269 (95% CrI 


state 



0-1134-0) 



Persistent 

//jp 

0-0252 (95% CrI 





0-1289-0) 



Transient 

Mu 

0-0684 (95% CrI 





0-266-0-0033) 



‘All’ represents infections without grouping of persistent and transient infections. Subscripts p, t, u and i represent persistent 
infection, transient infection, non-infectious state, and infectious state, respectively. CrI represents credible interval. 


infectious period, which is a reciprocal of recovery 
rate (y) (Table 2). Failure to detect coronavirus 
RNA in faeces or anal swabs may not necessarily 
imply recovery from the infection, and could conceiv¬ 
ably represent intermittent viral excretion or false¬ 
negative laboratory results. Flowever, recovery from 
SARS coronavirus infection in humans has been 
shown to occur when virus is no longer detected in 
faecal samples [30], In the absence of specific informa¬ 
tion about recovery from coronavirus infection in 
these bats, we assumed that failure to detect corona¬ 
virus in faecal samples or anal swabs similarly corre¬ 
sponded to recovery. 

The transmission rate (ft) could not be directly cal¬ 
culated from the CMR data analyses, and was instead 
calculated from basic reproduction number ( R 0 ) equa¬ 
tions. From Drexler et al.’s study [14] and from obser¬ 
vation of the CMR data [21], Smith [21] hypothesised 
that the initial epidemic peak was caused by the for¬ 
mation of a maternity roost of M. macropus from 
weeks 1 to 6, and that a second epidemic peak, after 
parturition, was caused by newborn pups who lost 
their passive immunity from weeks 6 to 12. Thus, 
with an assumption that a new epidemic began from 
week 7, we used an equation of R 0 - 1 + AD [31], 


where A represents the growth rate in an epidemic 
and D represents the average duration of the infectious 
period. We calculated A as 0-3328 per week by using 
the CMR data from weeks 7 to 12 (Supplementary 
figure SI). We estimated mean D as 1-7737 weeks, 
from the infectious period of recaptured bats (1/y) 
(Table 2). Thus, we estimated R 0 to be 1-5903. We 
estimated fi using R 0 - /?ND, where N is the total 
population size [31]. By assuming that N was 86 (the 
mean estimated size of the study population [21]) we 
could calculate /?, 0-0104. A rate of waning immunity 
(co) was unable to be estimated from this CMR data 
and has not been estimated in other studies. In the 
absence of data to suggest otherwise, we assumed 
the rate of waning immunity was comparable to that 
in human SARS-CoV infections [30]. Uncertainty in 
parameter values was included by sampling mortality 
and recovery rates from PERT distributions using the 
R package ‘mc2d’ [32], 

Epidemic model framework 

We built a deterministic density-dependent susceptible- 
infectious-recovered-susceptible (SIRS) model, using 
ODE (ordinary differential equations). We assumed 
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Fig. 2. Flow diagrams of epidemic models, (a) In the one-group model, three states were included: susceptible (S'), 
infectious ( I) and immune ( R ). (b) In two-group model, four states were included: susceptible ( S ), persistently infectious 
(Ip), transiently infectious (/,) and immune ( R ). Unlike the one-group model that included only one infectious state, the 
two-group model included two infectious states, and bats go through either a persistently infectious state or a transiently 
infectious state. The two infectious states have different recovery rates, in which the recovery rate (y p ) of a persistently 
infectious state is lower than the recovery rate (y t ) of a transiently infectious state, fr. transmission rate, y: recovery rate, co\ 
waning rate of immunity, /<: mortality rate, /: proportion of recaptured bats with persistent infection to the recaptured bats 
with persistent or transient infection. Subscripts p, t, u and i denote ‘persistent infection’, ‘transient infection’, 
‘uninfectious state’ and ‘infectious state’, respectively. 


that density dependence was appropriate for a corona- 
virus because SARS-CoV is transmitted in bats via a 
faecal-oral route, which is more suggestive of density- 
dependent transmission than frequency-dependent 
transmission [33, 34] and that SARS infection in 
humans has been modelled with density-dependent 
transmission [35, 36], Further, the number of bats in 
the roost was relatively small and this would allow 
homogenously mixed contacts among bats [21]. We 
used a SIRS model for coronavirus infection in bats, fol¬ 
lowing previous authors, who have used SIRS models 
for coronavirus infection in Miniopterus spp. [21]. 
Additionally, high SARS-CoV seroprevalence in bats 
supports the existence of a recovered class [4]. In the 
CMR data, positive detections of coronavirus RNA 
occasionally encompassed negative detections and 
negative detections occasionally encompassed positive 
detections, implying that infectiousness and non- 
infectiousness were not permanent [21]. 

To test effects of grouping of bats into persistently 
infectious and transiently infectious bats on the prob¬ 
ability of viral persistence in the bat population, we 
designed two alternative models: a ‘one-group’ model 
and a ‘two-group’ model (Fig. 2). The one-group 
model did not differentiate between persistently infec¬ 
tious and transiently infectious states, instead assumed 
that all recaptured bats (n = 23) with multiple or single 
detection of coronavirus RNA in the CMR data had the 
same recovery and mortality rates, regardless of the 
number of coronavirus RNA detections. On the other 
hand, the two-group model split bats into persistently 
infectious and transiently infectious groups. We 


assumed that recaptured bats with multiple detection 
of coronavirus RNA were ‘persistently infectious’ 
(77 — 7) and recaptured bats with single detection of cor¬ 
onavirus RNA in the CMR data were ‘transiently infec¬ 
tious’ (n - 16). Given the multistate model selection 
preferred the inclusion of the multistate effect 
(Table 3), we estimated different mortality rates for 
infectious bats (I class) and non-infectious bats (S and 
R classes). 

In the one-group model, the mean prevalence (P = 
0-2786) of coronavirus in the CMR data was used to 
set the initial number of infectious bats in the model 
with the population size (N = 86) and the remaining 
bats were considered susceptible (S-(l-P)N, I - 
PN, R = 0). In the two-group model, the initial num¬ 
ber of infectious bats was set based on the proportion 
of persistently infectious bats in infectious bats out of 
recaptured bats in the CMR data (f = 7/23) (.S’ = ( I P) 
N, I p = fPN, /, = (I f)PN and R = 0). 

Scenarios 

Six scenarios were set up based on different infectious 
periods. Scenario 1 was developed to describe the 
dynamics of coronavirus infection in M. macropus 
without grouping bats based on infectious period (the 
‘one-group’ model). Scenarios 2-6 were ‘two-group’ 
models. Scenario 2 was developed to describe the 
dynamics when bats were grouped into persistently 
infectious and transiently infectious bats, and scenarios 
3-6 were modifications of scenario 2, with extended per¬ 
iods of persistent infection. While we used the infectious 


Downloaded from https:/www.cambridge.org/core. Columbia University Libraries, on 30 May 2017 at 07:31:50, subject to the Cambridge Core terms of use, available at 

https:/www.cambridge.org/core/terms. https://doi.Org/10.1017/S0950268817000991 











6 J. Jeong and others 


Table 3. Multistate model selection with survival, recapture and transition probabilities of Myotis myotis 


Model 

number 

Survival probability 

Recapture probability 

Transition probability 

DIC* 

ADIC 1 

Parameter 

number 

Deviance 

1 

Group, multistate 

Group 

Group, multistate 

273-9 

0 

10 

258-2 

2 

Group 

Group 

Group, multistate 

274-1 

0-2 

8 

259-3 

3 

Multistate 


Multistate 

279-5 

5-6 

5 

267-3 

4 



Multistate 

280-1 

6-2 

4 

268-9 

5 

Group, multistate 

Group, multistate 

Group, multistate 

305-5 

31-6 

12 

251-1 

6 

Group 

Group, multistate 

Group, multistate 

319 

45-1 

10 

252-7 

7 

Multistate 

Multistate 

Multistate 

342-1 

68-2 

6 

258-4 

8 


Multistate 

Multistate 

352-7 

78-8 

5 

259-8 


CMR, capture-mark-recapture; DIC, deviance information criterion. 

Transition probabilities are probabilities that bats transit between infectious state and non-infectious state. Group means 
grouping of infectious bats into persistently infectious bats and transiently infectious bats. Multistate means the multistate 
effect of infectious and non-infectious states. Model 1 was found to be the most parsimonious model that best fits the 
CMR data. 

* DIC (a Bayesian analogy of the AIC (Akaike’s Information Criterion)) is a measure of the relative quality of statistical mod¬ 
els, and the model with the smallest DIC is estimated to be the model that would best fit the data. 

| ADIC is the change in DIC from the top-ranked model to each model. 


periods calculated from the CMR data analyses in 
scenarios 1 and 2, we assumed extended infectious 
period of persistently infectious bats to 5, 7, 9 and 
11 weeks in scenario 3-6, respectively, following a pre¬ 
vious study in which M. macropus could be identified 
with a putative novel Alphacoronavirus infection for 
up to 11 weeks [21], 

We simulated the model with 10 000 iterations, 
sampling from the range of parameter values calcu¬ 
lated from the CMR data and estimated the probability 
of viral persistence in the population of M. macropus. 
The R package ‘deSolve’ [37] was used to build the 
epidemic model. Time steps were weekly (following 
the time interval of the CMR data [21]) (see 
Supplementary material, R code S3). We assumed 
that the virus persisted in the population when at least 
one infectious bat remained at week 12. 

RESULTS 
CMR data analyses 

The multistate models enabled estimation of survival 
and recapture rates, of M. macropus, and transition 
rates between infectious and non-infectious states 
(Table 2). The survival rates were relatively constant 
during the period, whereas the recapture rates, gener¬ 
ally decreased in the first half and increased in the 
second half of capturing occasions (Fig. la). As a 
result, the transition rates should be considered with 
caution when recapture rate was low (for example, 
only one bat was captured at the fifth recapturing 


occasion) (Fig. lb). Paucity of the CMR data resulted 
in wide error bars. 

Multistate model selection 

The most parsimonious model supported grouping of 
bats into two groups based on the detection frequency 
of coronavirus RNA (single and multiple detections) 
for all three rates of survival, recapture and transition. 
There was little support for multistate effects on the 
survival rate (ADIC = 0-2, comparing models 1 and 
2; Table 3). This was in accordance with previous 
findings that coronavirus infection in bats showed no 
signs of illness [7, 8]. The results of model selection 
processes also indicated that the multistate effect 
should be excluded for recapture rate, and should be 
intrinsically considered for transition rate. 

Epidemic model simulations 

The SIRS epidemic models were simulated to generate 
the probabilities of coronavirus persistence in a bat 
population in the six scenarios. The simulated prob¬ 
ability of viral persistence in scenario 2 (two-group 
model, 0-5210) was somewhat higher than the prob¬ 
ability in scenario 1 (one-group model, 0-4094) 
(Table 4). This result showed that viruses were more 
likely to be maintained in the bat population when 
bats could either be persistently infectious or transi¬ 
ently infectious than when bats were not split into 
these groups. As the period of persistent infection 


Downloaded from https:/www.cambridge.org/core. Columbia University Libraries, on 30 May 2017 at 07:31:50, subject to the Cambridge Core terms of use, available at 

https:/www.cambridge.org/core/terms. https://doi.Org/10.1017/S0950268817000991 







Persistent infection of coronavirus 7 


Table 4. Probability of viral persistence based on varying periods of infection in one and two group models in six 
scenarios 



One-group model 

Two-group model 




Scenario number 

1 

2 

3 

4 

5 

6 

Transiently infectious period in weeks 
Persistently infectious period in weeks 
Probability of viral persistence 

1-774 (1/y)* 

0-430 

1- 165 (l/y t ) 

2- 982 (l/y p ) 
0-508 

1-165 (1/M 

5 

0-999 

1-165 (1/M 

7 

1 

M 65 (1/M 

9 

1 

1-165 (1/M 
11 

1 


The one-group model (with a single infectious period) was only used in scenario 1, while the two-group model (with both 
transient and persistent infectious periods) was used in scenarios from 2 to 6. Scenarios differed in the infectious period esti¬ 
mate used: scenario 1 and 2 used estimates from the CMR analyses (reciprocal of recovery rates (y), which are shown in 
Table 2), whereas scenarios 3-6 explored a series of hypothetical persistently infectious periods. The infectious period for tran¬ 
sient infections was constant for all two-group scenarios (based on 1/M, but the infectious period for persistent infections was 
varied. Specifically, scenario 2 used l/y p from the CMR analyses for persistently infectious period, whereas in scenarios 3-6, 
we assumed 5, 7, 9 and 11 weeks of persistently infectious periods respectively. 

* Scenario 1 used the same recovery rate for all bats. See Table 2. 


was lengthened to 5 weeks in scenario 3, the probabil¬ 
ity of persistence within the population reached to 
almost 1. In scenario 4-6, in which the periods of per¬ 
sistent infection were lengthened to 7, 9 and 11 weeks, 
respectively, the virus persisted in the population in 
every simulation. Although Smith [21] found that 
M. macropus can be identified with a putative novel 
Alphacoronavirus infection over periods of up to 11 
weeks, extension of the period of persistent infection 
from about 3 weeks (scenario 2) to 5 weeks (scenario 3) 
markedly increased the probability of viral maintenance 
to almost one. 

DISCUSSION 

This study analysed CMR data of a species of 
Australian bats (M. macropus) with a putative novel 
Alphacoronavirus infection in order to decide whether 
it is appropriate to divide bats into persistently infec¬ 
tious bats and transiently infectious bats, and to deter¬ 
mine whether viral maintenance is improved by the 
differentiation of bats into these groups. Grouping 
of bats into persistently infectious bats and transiently 
infectious bats was supported by the CMR multistate 
model selection processes (Table 3). By exploring vari¬ 
ous scenarios in epidemic models, we found that 
population-level viral persistence was marginally more 
probable when infectious periods were heterogeneous 
(bats were either transiently or persistently infectious; 
Table 4). In addition to a previous study that suggested 
persistent infection of coronavirus in bats in North 
America through the small number of bats sampled 
[23], our study provides much validity by using larger 


number of bats sampled [21], although the sample size 
is still not large enough to provide robust estimates. 

Further understanding of the potential and likeli¬ 
hood of within-host persistent infections is important 
to understand spillover of bat-borne viruses. Persistent 
infections contribute not only in maintaining viruses 
in a population but also in triggering pulses of virus 
shedding. Temporally synchronised stressors such as 
food shortage and reproduction may weaken the 
immune system of bats, facilitating persistent infec¬ 
tions in infectious bats [38], Thus, those temporally 
synchronised stressors can cause viral shedding from 
persistently infectious bats during a limited period, 
leading to pulses of viral shedding. Because spillover 
of bat-borne viruses has been associated with pulses 
of viral shedding from bats [38], the finding of persist¬ 
ent infection in bats contributes in understanding spill¬ 
over mechanisms of bat-borne viruses. 

Maternity roosts were assumed to play an import¬ 
ant role in maintaining coronavirus at the population 
level [24]. In coronavirus transmission within a popu¬ 
lation of Myotis, an epidemic peak was observed at 
the formation of the maternity colony and another 
peak was observed as newborn pups lost their passive 
immunity [14], Transmission of viruses during the 
breeding season is more readily facilitated, compared 
with other times due to frequent contact among indivi¬ 
duals within a maternity roost [14, 21, 39], Fluctuating 
recapture rates and coronavirus prevalence are pre¬ 
sumed to be related to the changing ethology of bats 
during the breeding season. Therefore, CMR data cov¬ 
ering both breeding seasons and non-breeding seasons 
are required to investigate how coronaviruses can be 
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maintained in a bat population even if maternal roosts 
do not facilitate viral transmission. 

These findings need to be treated with some cau¬ 
tion, because the low sample size limits our confidence 
in the results of modelling. For example, the recapture 
rates were markedly different from recapturing 
occasion to occasion (Fig. la), as only one bat was 
sampled at the fifth recapturing occasion [21]. 
Additionally, the CMR data included only corona- 
virus RNA detections, without detections of corona- 
virus antibodies. A SARS coronavirus crude antigen 
ELISA was not successful in detecting antibodies 
against the putative novel Alphacoronavirus in these 
bats [21]. Thus, the lack of seroprevalence in the CMR 
data limited the accuracy of the estimated duration of 
recovery. Nevertheless, the CMR data deserve intensive 
analyses because longitudinal sampling of individual 
bats for infectious diseases is difficult to achieve and 
allows estimation of epidemiological parameters. The 
CMR study in Smith's study [21] had relatively high 
recapture rates thanks to finding a lifting hole of a bridge 
where bats, infected with coronavirus, roosted with a 
high affinity to their colony. Fluctuating recapture 
rates across capturing occasions in the CMR data, 
which might be associated with trap-shy and adaptive 
methods of capturing, highlighted the difficulty of recap¬ 
turing specific bats during the entire sampling period. 

Persistent infections have been hypothesised as a 
mechanism for viral persistence, not only in corona- 
viruses [23], but also in other bat-borne viruses, for 
example filoviruses and henipaviruses, despite lack 
of direct evidence [20], It was also hypothesised that 
a small portion of super-long-shedder bats with a 
long infectious period could maintain Hendra virus 
in a bat population [13], Although further comprehen¬ 
sive datasets are needed to understand the effects of 
persistent infection, this study shows the value of longitu¬ 
dinal individual data from bats with viral infections to 
elucidate transmission dynamics of bat-bome viruses, 
underscoring the need of individual bat tracking data 
with infection states to improve our understanding of 
infection dynamics of bat-bome viruses. 

SUPPLEMENTARY MATERIAL 

The supplementary material for this article can be 
found at https://doi.org/10.1017/S0950268817000991 
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